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Abstract 

The kinetic instabilities of the Weibel-type are presently invoked in a large variety of astrophysical 
scenarios because anisotropic plasma structures are ubiquitous in space. The Weibel instability is 
driven by a temperature anisotropy which is commonly modeled by a bi-axis distribution function, 
such as a bi-Maxwellian or a generalized bi-Kappa. Previous studies have been limited to a bi- 
Kappa distribution and found a suppression of this instability in the presence of suprathermal tails. 
In the present paper it is shown that the Weibel growth rate is rather more sensitive to the shape 
of the anisotropic distribution function. In order to illustrate the distinguishing properties of this 
instability a product-hi- Kappa distribution is introduced, with the advantage that this distribution 
function enables the use of different values of the spectral index in the two directions, Ky 7^ n±. 
The growth rates and the instability threshold are derived and contrasted with those for a simple 
bi-Kappa and a bi-Maxwellian. Thus, while the maximum growth rates reached at the saturation 
are found to be higher, the threshold is drastically reduced making the anisotropic product-bi- 
Kappa (with small kappas) highly susceptible to the Weibel instability. This effect could also 
rise questions on the temperature or the temperature anisotropy that seems to be not an exclusive 
source of free energy for this instability, and definition of these notions for such Kappa distributions 
must probably be reconsidered. 
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I. INTRODUCTION 



The huge amount of incoming radiation and ionized particles from space suggests the 
wide-spread existence of kinetic anisotropics in cosmic plasmas. Moreover, space plasmas 
are sufficiently dilute and their collisionality is sufficiently low (see, e.g., table 8.1 in 
and, therefore, expected to be far from a Maxwellian equilibrium (only provided by the 
short-range binary collisions of plasma particles). Thus, we have direct proofs from the 
observations and in-situ measurements that suprathermal populations are widely present 
at different altitudes in the solar wind plasma [2 and probably in the solar corona jo]. 
These populations exhibit suprathermal high energy tails and are fitted quite well by the 
family of Kappa distribution functions [3], which are power laws in particle speed. The 
Kappa distribution generalizes the notion of equilibrium in coUisionless plasmas deviated 
from a thermal (Maxwellian) equilibrium, and which contain fully developed turbulence in a 

on the basis of a superadditive entropy [9|, 110| . 

The so-called Kappa distributions are defined by using a spectral index k, which deter- 
mines the slope of the high energy tails in the velocity spectrum of plasma particles. In 
the limit of very large k — > oo, the Kappa functions degenerate into Maxwellians IjJ. The 
effects of these suprathermal populations and their anisotropics on the threshold conditions 
and the linear growth rates of kinetic instabilities have been studied by modeling plasma 



with diverse anisotropic distributions fully or partially populated by Kappa particles 



12 



27l |. These effects largely vary, depending on the shape of the distribution function and the 



nature and the frequency of the plasma mode. A selection of Kappa distributions and their 
limiting forms for a large k — )■ oo is given in Ref. [ll|. Table I. The most simple anisotropic 
distribution is the bi-Kappa distribution function introduced here in Eq. ([1]), and which 
describes a temperature anisotropy, T± ^ Ty, with the same index k\\ = = k for both 
directions. This distribution has widely been used for describing kinetic instabilities in a 
drifting or a non-drifting plasma, e.g, the ion acoustic instability 12|, the electromagnetic 



ion-cyclotron and firehose instabilities [1 



and electron eye 



Weibel-type 
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otron instabilities 
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Recently, Basu 




, the mirror unstable mode 16|, the whistler 



24| . and the nonresonant instabilities of 



27|. 



26| has reviewed the stability properties of hydromagnetic waves in a 
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plasma distributed after a product-bi-Kappa function such as the one introduced here in 
Eq. dlj). The product-bi-Kappa distribution function represents a generahzation of the 
Kappa-Maxwelhan distribution, which is a product of a one-dimensional Kappa distribution 
along a preferred direction in space, e.g^ the guiding magnetic field, and a Maxwellian 
distribution in the perpendicular plane 17|. In a hybrid Kappa-Maxwellian plasma, unlike 
a uniform Maxwellian or a uniform Kappa, the dispersion properties and the stability were 



found to be markedly changed 



11, 



19 



201 ] ■ Notable is that, because of the anisotropy of the 



contours in the velocity space, such a Kappa-Maxwellian distribution can be unstable even 
for equal parallel and perpendicular temperatures, or some very well-known unstable modes 
become stable to a temperature anisotropy in such asymmetric distributions 20] . 

The product-bi-Kappa distribution function has not often been applied so far, but we 
think that this function deserves further consideration for the following reasons: (1) it 
decouples the dynamics of the plasma particles over the two principal directions allowing 
not only for distinct temperatures, Ty and Tj_, but for distinct spectral indices, k\\ and as 
well; (2) there is a possibility to regulate this coupling in accord to a guiding magnetic field 
or other external constraints by the interplay of k,\\ and k±; (3) compared to a bi-Kappa, 
a product-bi-Kappa distribution allows for a more realistic description and (4) enlarges the 
number of distributions with an elaborated dispersion approach; (5) one more important 
feature (recently revealed for a similar Kappa-Maxwellian plasma 20|) is that such a product- 
bi-Kappa model permits further analytical progress leading to tractable expressions for the 
dielectric tensor elements enabling, for example, the investigation of oblique modes in a 
magnetized plasma. 

Electromagnetic instabilities of the Weibel-type 
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29| are driven by an arbitrary devia- 



tion of the particle velocity distribution from its equilibrium, whether it is a bulk (relative) 
motion of streaming particles 29|] or a temperature anisotropy 28| . Kinetic anisotropics are 
ubiquitous in space plasmas (extending from heating flows and temperature anisotropics to 
particle jets and shock waves in outflows, or interpenetrating plasma shells in interplane- 
tary wind) and, presently, Weibel instabilities are invoked in a large variety of astrophysical 
scenarios. These instabilities create magnetic field such as the cosmological seeds necessary 
for the dynamo mechanism 30|, l3l| , or a quasi-stationary magnetic field boost for the syn- 



chrotron emissions in astrophysical sources 



32| . Moreover, the magnetic field fluctuations 



observed in interplanetary space are also attributed to these instabilities (e.g., filamenta- 
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tion, whistler, mirror, oblique firehose) 33|, |3J], which compete with other constraints of 



plasma particles (adiabatic expansion. Coulomb collisions) and maintain a relatively small 



temperature anisotropy in the solar wind 



34 



38 



Here we discuss the Weibel instability driven by a temperature anisotropy of plasma 
particles. Recent studies of this instability have been limited to a bi-Kappa distribution and 



found that the instability is suppressed in the presence of the Kappa tails [2l|, |23| . For the 
same arguments formulated above, we investigate here the product-bi-Kappa distribution as 
a favorable alternative that seems to add an excess of free energy in the velocity space, and 
to enhance the instability. In the present paper we proceed to a comparative analysis of the 
effects of these two distributions, the bi-Kappa and the product-bi-Kappa, on the Weibel 
instability. The plasma is assumed to be coUisionless and spatially homogeneous. 

In order to describe the initially unperturbed plasma system we first introduce the bi- 
Kappa distribution function 

using polar coordinates {vx,Vy,Vz) = (f± cos 0, sin 0, f y) in the particle velocity space. 
This distribution is normalized to unity, J d^v Fi = 1, and makes reference to the equivalent 
thermal velocities 0\\,±, which relate to the effective temperatures of the plasma particles: 

^11 - ^ = f / M^i(-lh-x) = ^ i^Of (2) 
Ikb Kb J 2kB 2k — S " 



2kB 2kBj " 2kB2K-2> ^' ^ ' 

for a spectral index n > 3/2. 

The second distribution function examined here is a product-bi-Kappa 



2 \ -''11-1 / 2 \ -Kx-1 



as was defined in Summers and Thorne This function is also normalized to unity 

J d^v F2 = 1, and the equivalent thermal velocities 6\\^± are given by 
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FIG. 1: Contours of the distribution functions ([T|) and dH are plotted with dotted and dashed lines, 
respectively, in panel (a) for = kj_ = n = 3 and vt^/c = 2vt^Jc = 0.2, and in panel (b) for a 
very large k — )• oo leading to an exact fit of these two contours with that for a bi-Maxwellian (solid 
line). A Kappa-Maxwellian distribution function Q with = 3 and ^ oo (dot-dashed line) 
is plotted in (c), and all these distributions are compared in panel (d). For a small Kappa, remark 
the high energy tails for both distributions and the prominent anisotropy of the distribution (jl]). 

for different spectral indices k,± > 1 and k,\\ > 1/2, respectively. 

Contour plots of the distribution functions ([T]) and are displayed in Fig. [H At low 
values of k, there are significant differences between these two distributions and a prominent 
asymmetry of the product-bi-Kappa distribution, which indicate a surplus of temperature 
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anisotropy, e.g., in panels (a) and (d). However, both these distributions functions approach 
the same bi-Maxwelhan in the hmit of a very large spectral index, e.g., in panel (b). 

Assuming an excess of perpendicular temperature (T± > Ty), the Weibel instability 
develops along the parallel direction (subscript " ||") with a wave-number k = k^\, aperiodic, 
Re(c<;)= cOr = 0, and with a growth rate Im((X')= > . By using standard techniques 
based on the linearized Vlasov-Maxwell equations, it is easy to find the dispersion relation 



for the electromagnetic modes propagating along the parallel direction 



21 



23| 



k^c^ 



l + Ttk 



dv\\ 



oj — kvw 



^ 3 dF ^ 



(7) 



where ojp = {Airnee^ /nieY^'^ is the plasma frequency. 

First we insert the bi-Kappa distribution function ([1]) in Eq. ([7]) and find a dispersion 



relation 
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23| 



u — k^c 



2^2 



oj — k^c 



2^2 



1 + 


01 

el 


1 + 


kO\\ 




ykoi 


1 + 


7]| 


1 + 


UJ 

kOi\ 




i- 



(8) 



in terms of the modified plasma dispersion function 

1 m 



53(/) > 0. 



(9) 



7ri/2fi;i/2 T[k-1 /2] x-f 
This dispersion relation admits purely growing solutions of the Weibel type for wave-numbers 



smaller than a cutoff value 



1/2 



UJr, 



Til 



- 1 



1/2 



(10) 



The exact numerical growth rates are displayed in Fig. [2] with solid lines for small values of 
K, and with dotted lines for a very large k — ?■ oo (bi-Maxwellian plasma). In this case the 
cutoff wave- number does not depend on the spectral index k > 3/2 and this is confirmed 

1 /2 

in Fig. [2] where kdc/up = (Tj_/T|| — l) ^ 1.73 is the same for all three cases. For the 
existence of a finite kd 7^ 0, the instability threshold r = min(ri/r||) is simply found as 
T^/Tii >r = l. 

If we insert the second distribution function (j4j), the dispersion relation ((Tj) takes a new 
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FIG. 2: The growth rates of the Weibel instabiUty, solutions of Eqs. ([8]) (sohd hnes) and (jlip 
(dashed hnes) for an anisotropic Kappa distributed plasma with vt^/vt^^ = 2, vt^^ = 0.02c and 
for three values of K|| = k = 2 (in a), k =4 (in b), and k = 10 (in c). The growth rates for a 
bi-Maxwellian plasma {k — )■ oo) are plotted with dotted lines. 



and different form 
= 



- k^c^ ^ el k^{k\\ + 1/2) 
0] /t||(/t±-l) 



UJt. 



- /til + 1/2 
1 + ^^^ 



where 



z.Af) 



1 



T|| Kx. - 1/2 

rhi + 2] 



1 + — , 

kOw " V ke 



(11) 



(12) 



vtVSkJ/^ r + 3/2] X - / 

This dispersion relation does not involve any dependence on but only on Ky, and it 
can therefore be further applied to describe the instability of a Kappa-Maxwellian (one- 
dimensional Kappa distribution in the parallel direction and a Maxwellian distribution in 



7 



the perpendicular plane). Notice that, both the plasma dispersion functions ([9]) and ( !T2|) 
approach the standard dispersion function of Fried and Conte jsol in the limit of a very large 

K, K|| — )• OO. 

The aperiodic solutions of dispersion relation f ITT]) are displayed in Fig. |2] with dashed 
lines. In this case not only the growth rates increase (as was expected from Fig. [H where the 
excess of anisotropy of a product-bi-Kappa distribution compared to a bi-Kappa, is evident), 
but the instability extends as well to larger wave-numbers (smaller wave-lengths) due to a 
non-negligible dependence on the spectral index k\\ of the new cutoff wave-number 

1/2 



k 



c2 



Ti_ Kii + 1/2 



1 



(13) 



c [T|| K|| - 1/2 

If we look, for example, at the growth rates displayed with dashed line in Fig. |2^, and 
which are solutions of Eq. f lTTj) for an index k\\ = k = 2, and a temperature anisotropy 
vt^/v^ = 2 (where f ^| = -y/2/cBTij|/m), the cutoff wave-number is indeed larger kc2c/ujp = 



54, /{3n 



1/2 



2.38 > kci — 1.73. The cutoff wave-number ( !T3l) decreases and 



approaches ( ITOl) only for a sufficiently large k — > 00. 



Unstable bi-Kappa distributions 



Unstable 




FIG. 3: Contrast of the instability thresholds showing the unstable bi-Kappa distribution in the 
gray (orange online) region in the limits of > 1.5 and T±/T^^ > 1, and the unstable product- 
bi-Kappa distribution in the light-gray (light-orange online) region limited by > 0.5 and the 
condition (|14p for the temperature anisotropy, and including the more restrained gray region. 



Furthermore, for a product-bi-Kappa, the anisotropy threshold r decreases to less than 
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unity 




K|| - 1/2 



< 1. 



(14) 



K|| + 1/2 



and towards the limit of k,\\ — i- 1/2, the threshold vanishes, r ~ 0, and the instability can 
grow freely. If = 2 the instability will develop for T±/T\\ > t = 3/5, and then the 
instability threshold increases asymptotically to unity, r — )■ 1, as the spectral index becomes 
very large k\\ — )■ oo. 

To conclude, in this paper we have revised the effects of Kappa anisotropic distributed 
plasmas on the Weibel instability. Previous models were limited to a bi-Kappa distribution 
and have been extended by introducing a more general product-bi- Kappa distribution. In 
this new distribution, the dynamics of plasma particles over the two principal directions are 
decoupled and characterized by two distinct temperatures and two distinct spectral indices, 
K,± 7^ K,\\. The growth rates and the instability threshold have been found to be very sensitive 
to the shape of the anisotropic distribution. While for a bi-Kappa distribution the Weibel 
instability is suppressed, for a product-bi-Kappa the growth rates are enhanced and the 
instability threshold is significantly lowered by comparison to a Maxwellian. 

The enhancing effect can not be attributed to a resonant interaction with the energetic 
particles from Kappa tails because this instability is nonresonant, but it is fully supported by 
an excess of anisotropy and free energy of the product-bi-Kappa distribution in the velocity 
space. In this sense, the contour plots from Fig. [1] are very suggestive: if we compare 
contours in panels (a) and (b), the bi-Kappa is less asymmetric than the bi-Maxwellian 
leading to lower growth rates of the Weibel instability, while the product-bi-Kappa is more 
asymmetric than the bi-Maxwellian and enhances the instability. 

The temperature anisotropy threshold for the Weibel instability can be drastically re- 
duced in anisotropic plasmas with product-bi-Kappa distribution functions that suggests 
two possible explanations: either the fundamental notions, such as the temperature or the 
temperature anisotropy must be redefined for these Kappa distributions, or the temperature 
anisotropy is not an exclusive source of free energy for this instability. 

These results can easily be extended to high beta plasmas (/3 = thermal energy/magnetic 
energy) widely present in space and where the ambient magnetic field has only a minor 
influence without changing the essential features of the instability, but sustaining such 
asymmetric distributions with k± > due to some equilibration and isotropization in the 
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perpendicular plane, and a preferential motion and acceleration along the magnetic field. 
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